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Abstract 

The two-phase composite approach of Estrin et al. (1998) describes an 
evolving dislocation cell structure. Mckenzie et al. (2007) enhanced the model 
to capture the effects of hydrostatic pressure and temperature during severe 
plastic deformation. The goal of the present study is to incorporate this 
microstructural model into the macroscopic viscoplasticity framework pro- 
posed by Shutov and KreiBig (2008a). Thereby, the two-phase composite 
approach is examined carefully. Both physical and numerical drawbacks are 
revealed and possible solutions are presented, thus leading to a refined micro 
model. Moreover, some improvements concerning reliable parameter identi- 
fication are suggested as well. The material parameters of the refined micro 
model are identified for an aluminum alloy using TEM cell size measure- 
ments. Then, an extension to non-proportional deformation is performed 
in such a way that the evolution of dislocation densities becomes sensitive 
to load path changes. Experimental findings suggest that such deformation 
modes can significantly influence the evolution of microstructure, including 
the dissolution of cells and the reduction of total dislocation density shortly 
after the load path change. In order to capture such effects, some tensor- 
valued state variables are introduced which couple the refined micro model 
with the macroscopic viscoplasticity model. As a result, a new system of 
constitutive equations is obtained. In order to demonstrate its capability to 
respond to load path changes, load CclSGS cLS typical for Equal Channel Angu- 
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lar Pressing (ECAP) are considered. The obtained evolution of dislocation 
populations differs signficantly depending on which ECAP route is applied. 

Keywords: Two-phase composite model, dislocation cells, viscoplasticity, 
non-proportional loading, ECAP 



1. Introduction 

An adequate description of evolving dislocation structures is of great im- 
portance for physically based modeling of many metallic materials. There 
are sophisticated multi-scale approaches to model dislocation dynamics in a 
discrete or continuous manner (e.g. Devincre et al., 2001; Groh et al., 2009; 
Cai et al., 2006; Le and Sembiring, 2008). Thus, precious insights into the 
collective behaviour of dislocations can be obtained (for an overwiew, see Bu- 
latov and Cai, 2006). However, the scale bridging to engineering plasticity 
applications where large deformations and long time frames occur is still an 
open problem. 

The number of degrees of freedom can be drastically decreased by intro- 
ducing dislocation densities as state variables^ into continuum scale models. 
Corresponding evolution equations can be deduced from microstructural con- 
siderations. Depending on the number of introduced dislocation populations 
these approaches are called one-parameter- mo del (only total dislocation den- 
sity) or multi-parameter-models, respectivly (cf. Goerdeler and Gottstein, 
2001; Ding and Guo, 2001). Within the last decade, this approach has been 
used by many researchers to include an evolving dislocation structure in 
continuum- scale modeling of metal plasticity. In particular, such informa- 
tion can be incorporated within models of crystal plasticity (see Gerard et al., 
2013, and references therein) or even purely phenomenological modeling (see 
Shutov and Thiemann, 2012). 

There are numerous theoretical and experimental indications that dis- 
locations within crystals constitute self-organizing systems (Amodeo and 
Ghonicm, 1990b; Yamaki et al., 2006; Grcgor, 1998). As a characteristic fea- 
ture, they tend to distribute inhomogeneously and thus form spatial patterns. 
Among others, the three-dimensional cell structure (Langlois and Berveiller, 
2003) is typical for a class of metals and alloys. After the dislocation struc- 
ture has emerged, further monotonic deformation results in the shrinkage 
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of the structural characteristics: dislocation spacing, wall thickness and cell 
size. In the case of non-monotonic deformation the behavior seems to be 
different. Experimental findings (Hasegawa et al., 1975, 1986) suggest that 
load path changes can lead to the dissolution or disruption of cells. 

The evolving cell structure, i.e. the formation, sharpening and possible 
dissolution of the cells, governs the strength behaviour on the macroscopic 
level. Additionally, dislocation cells can serve as precursor of forthcoming 
subgrains even leading to nanocrystalline microstructures under severe plas- 
tic deformations (Toth et al., 2010). In particular. Equal Channel Angular 
Pressing (ECAP) is a promising technique which allows development of an 
ultrafine-grained microstructure without net change in the billet's shape (Se- 
gal, 1995). 

The present study is organized as follows: the two-phase composite ap- 
proach of Estrin et al. (1998) is outlined and carefully analyzed (cf. Section 
2). Both physical and numerical drawbacks are identified and possible solu- 
tions are presented, thus leading to a refined model. The problem of a rehable 
parameter identification is discussed in Section 3. There, material param- 
eters of the refined model are identified for the aluminum alloy AA 6016. 
In Section 4 the model is extended to non-proportional deformation in such 
a way that the evolution of dislocation densities becomes sensitive to load 
path changes. Towards that end, the refined micro model is coupled with 
the macroscopic viscoplasticity model of Shutov and Kreifiig (2008a), which 
is also outlined briefiy. Finally, a new system of constitutive equations is ob- 
tained including the micro model of Mckenzie et al. (2007) as a special case. 
In Section 5, the evolution of dislocation populations during Equal Channel 
Angular Pressing (ECAP) is considered. In this context, the extension to 
non- proportional loading is especially important, since most commonly used 
ECAP routes involve distinct load path changes. Finally, simulation results 
for routes A, C, Be, and E are discussed. 

Let us briefiy introduce the notations used in this paper. A second rank 
tensor is denoted by a small or capital letter with two underscores, e.g. U. 
The coefficients with respect to a certain Cartesian coordinate system 
are Uab = ■ U ■ Ch- As usual, the trace of a tensor is denoted by the 
operator tr. This enables the definition of the double contraction of two 
second-rank tensors hy U •• V — tr(C/ • V"). For the scalar product another 
double contraction in the form U : V — tr( [/ • V"""^) is used. Arranging 
second-rank tensor coefficients in a quadratic matrix is denoted by [Uab] or 
just [U]. Bold symbols indicate real- valued matrices regardless of where the 
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coefficients arise from, e.g. x = [1, 2, 3, 4]'^ is a 4-tuple or column matrix. 

2. Two-phase composite model 

Tfie fundamental assumption of the dislocation density based two-phase 
composite model of Estrin et al. (1998) and Mckenzie et al. (2007) is that a 
dislocation cell structure has already formed. Among the variety of reported 
dislocation configurations (Amodeo and Ghonicm, 1990a), cells are charac- 
terized by regions of lower dislocation density (the cell interior) surrounded 
by high dislocation density walls. The latter represent a topologically contin- 
uous skeleton (Mckenzie et al., 2007). Hence, the polycrystal can be thought 
of as a solid consisting of grains which again consist of dislocation cells. 
Because of the large difference in dislocation density, the cell wall and the 
interior show different mechanical behaviour. Thus, the solid can effectivly 
be treated as a two- phase composite (Mughrabi, 1987). Whereas the effect 
of grain orientations can be incorporated by using crystal plasticity (e.g. Yal- 
cinkaya et al., 2009), further modeling is necessary to take into account the 
two phases due to dislocation cells. 

2.1. Outline of the existing model 

The "unit cell" of the dislocation arrangement is a single cube with an 
edge length of d and a wall thickness of u)/2 ^. The cell wall volume over 
the total volume Vt yields the volume fraction / of cell walls: 

= ^ id-wr^d\i-f). (1) 

As a modeling assumption, the dislocations are classified with respect to their 
position. A part of the total dislocation length Lt is arranged within the 
cell wall volume and the remaining part is located in the cell interior 
Vc- Hence, we can define the densities of two dislocation populations: 

p^^— , p^ = — . (2) 



^In some studies, spherical dislocation cells are assumed in order to simplify the stress 
computation (Sauzay, 2008; Brahme et al., 2011). 
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Figure 1: Scheme of the dislocation ceh arrangement, the "unit ceU" is colored black. 

Dividing the total dislocation length Lt = + L^. by the total volume 
Vt — V^-\-Vc yields the rule of mixtures for the total dislocation density: 

Pt = /Pw + (1 - /)Pc . (3) 

By Eq. (3) it is not explicitly determined how the populations behave, e.g. 
whether they are mobile or immobile^. However, Estrin et al. (1998) consid- 
ered the dislocations woven into the walls as immobilized. This point will be 
important for Section 4.3. 

By now, no evolution of cell size d has been included explicitely. In the 
geometrical model (cf. Fig. 1) w represents a fraction of the cell size d. Thus, 
changing the wall thickness w or volume fraction / docs not affect the cell 
size. According to Amodeo and Ghoniem (1988) and Estrin et al. (1998), it 



■^There exist alternative criteria to group dislocations, e.g. by their character (edge or 
screw) or mobility (mobile or immobile) (e.g. Amodeo and Ghoniem, 1988; Austin and 
McDowell. 2011, 2012; Gao and Zhang, 2012). Thus, one rather yields chfferent species 
(dislocations with the same behavior) instead of populations (dislocations in the same 
domain). In contrast to Eq. (3), the total dislocation density is then just the sum of the 
densities of species and there is no rule of mixtures. 
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is assumed that the total dislocation density determines the cell size by the 
empirical relation 

d = K/^, , (4) 

with a proportionality constant K. In order to derive evolution equations 
for the dislocation densities in the cell wall and the interior, elementary 
dislocation interactions are considered. Finally, one obtains the following 
system of evolution equations (cf. Estrin et al., 1998; Mckenzie et al., 2007): 



dpw _ ^ 6/3*7c(l - ff'^ ^ \/3/3*7c(l - /) r- 



dt bdf bf 



- ko exp — - TwPw , (5a) 

- H 7^\fP^ - ^0 — 7cPc ■ (5b) 



dt K1-/)V3 VA^w .^U^^^ 

The terms on the right hand side represent the following processes: flux of 
dislocations from cell interior into the walls with ^* representing the fraction 
of Lc woven into the wall (1^* term in (5a) and (5b) respectivly) ; generation 
due to multiplication by Prank- Read sources (2°*^ term in (5a), (5b)) with 
a* as fraction of active sources; and annihilation supported by climb of edge 
dislocations (S""*^ terms in (5a), (5b)). All these processes'* are related to 
a change in and respectivly, which contributes to the change in 
and Pc- The shear rates in the cell wall and the interior are denoted by 7^ 
and 7c respectively. The remaining parameters are explained in Table 1. It 
can be argued that the "stiff" cell walls impose the same deformation on 
the cell interior (Estrin et al., 1998) leading to the kinematic assumption 
7w = 7c = 7- Thus, the evolution equations and the initial conditions can be 



Note that the nucleation of new dislocations is not considered within this model. 



6 



Table 1: Parameters of the dislocation density based two-phase composite model from 
Estrin et al. (1998) 



Para. 


Value/ Range 


Brief explanation 


h 


2.86-10-^°m 


magnitude of Burgers vector for Al 


Q 


1.2-10-29m3 


atomic volume for Al 




1.38065- 10-23J/K 


Boltzmann's constant 


a 


> 


controls the evolution of / (cf. Eq. (7), (49)) 


foi foo 


(0,1) 


initial and final values for / 


P 


> 


hydrostatic pressure, process dependent 


T 


> 


absolute temperature, process dependent 


7 


> 


plastic shear strain rate, process dependent 


7o 


> 


reference shear strain rate 


a* 


[0,1] 


fraction of active Prank-Read sources 




[0,1] 


emigrating fraction of pc 




> 


dynamic recovery exponents cell interior/wall 




> 


dynamic recovery constant 


K 


> 


proportionality between d and 1 / 



written as 



dp. 6/3-(l-/)^/3 ^ , \/3/3*(l-/) ^ 
~ ZTTl vPt H TT VPw 



d7 hKf • bf 



ko exp — - — Pw , (6a) 

n^k^T J V7o/ 



dpc 6/3* ^ 6a* ^ /7 



-l/nc 

Pw|7=0 = PwO , Pel 7=0 = PcO • (6c) 

Mathematically, (6a) and (6b) constitute a set of two coupled, nonlinear, 
ordinary, first order differential equations for the variables p. and pc- Since 
the wall volume fraction f{t) is time-dependent, system (6) is rheonome. 
Assuming the evolution of the cell wall volume fraction in the form (Mckenzie 
et al., 2007) 

/ = /oo + (/o - foo) exp(-Q;s) , (7) 

where s denotes the accumulated plastic strain, initial value problem (6) can 
be solved numerically. Consideration of the velocities [dpw/d7, dpc/d7]^ in 
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the phase space [p^,pc]^ for constant values of / reveals a fixed point char- 
acteristic. In Fig. 2 the fixed point is characterized by a vanishing velocity. 
In particular, system (6) exhibits fading memory behavior. Moreover, the 
solution is sensitive to the evolution of / (cf. Fig. 2 and 3). 
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Figure 2: Sketch of the velocity field in the phase space for / = /oo = 0.05 (remaining 
parameters taken from Tabic 1). The circles represent the grid. The zoom box on the 
right side reveals the fixed point. 



2.2. Critical analysis, improvements and modifications 

The model of Estrin et al. (1998) has been widely used in literature and 
adopted to many applications, especially in the field of severe plastic defor- 
mation (Baik, 2003; Hosseini and Kazeminezhad, 2009; Toth et al., 2010). It 
has to be emphasized that the derivation of the evolution equations is based 
on clear physical considerations of few elementary dislocation interactions. 
Still, the model predictions are applicable over a wide range of strain and fit 
experimental results very well (Mckenzie et al., 2007). 

Nevertheless, the model also has some weak points. It is known (Madec, 
2002) that cellular patterns form at the onset of stage III of the deformation 
of single crystals or at small plastic strains in multislip conditions. In any 
case, they do not emerge until the accumulated plastic strain reaches some 
value Scrit > (Chinh et al., 2004). Strictly speaking, the model is not valid 
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before since it cannot take into account the precedent self-organization of 
dislocation structures. At the early stage of plastic deformation, alternative 
models should be implemented which account for the influence of the grain 
boundary topology on the dislocation storage (cf. Delannay and Barnett, 
2012). Further, the characteristic microstructural quantities, such as the cell 
wall volume and the dislocation densities and are not clearly defined 
for s < Scrit since they do not exist^. Thus, it is unclear how to choose 
initial values for them. It is questionable whether measurements from the 
undeformed sample can be used. Fortunately, initial values are important 
only at the beginning due to the fading memory of the model. This aspect 
and a possible treatment are discussed in Section 5. 

Next, as seen in Fig. 2 and 3 the solution of system (6) is very sensitive 
to the cell wall volume fraction /. Consequently, the evolution of / should 
be measured accurately. However, this is related to huge experimental effort 
and not always justifiable in practise. 

A minor but noticible problem stems from the treatment of the constant 
7o. As it is merely used to produce a dimcnsionlcss term and set to 1/s in 
(Mckenzie et al., 2007) one should expect that the solution of system (6) is 
independent from its choice. Because 7/70 appears in (6a) and (6b), but with 
two different exponents n^, ric, the solution does depend on 70. Indeed, given 
a set of parameters ko,ny^, ric which correspond to a certain 70, no equivalent 
parameters can be found if 70 is replaced by some new value. Hence, it should 
be considered as an additional parameter to be identified. 

As an evident requirement, a physical model should produce reasonable 
results under any set of parameters within their range, but it was found that 
system (6) can lead to negative values of pc, e.g. using the parameters given 
from Mckenzie et al. (2007) and a — 0.1. This tendency is visible in Fig. 3 
where the velocity vectors point downwards in the vicinity of — ?■ 0. If a 
trajectory begins in these regions of phase space it is likely to cross pc = 
and go below. Obviously, the definition of dislocation density (cf. Eq. (2)) 
requires it to be non-negative. 

Therefore, one of the goals of the current publication is to modify the 
original evolution equations, such that the positiveness of dislocation densi- 



^Note that some authors suggest the initial condition / = 1 in order to describe the 
state when no dislocation pattern is present (Yalcinkaya et al., 2009). 
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10^' 10^' 10^' 

Figure 3: Sketch of the velocity field in the phase space for / = /o = 0.07 (remaining 
parameters taken from Table 1). 

ties will be guaranteed^. A closer look at the derivation from Estrin et al. 
(1998) reveals that the problem originates from the use of Orowan's relation. 
More precisely, this equation is applied to relate the mean velocity Vc of mo- 
bile dislocations in the cell interior with the macroscopic shear strain rate 
7c: 

^^c = 7c/ {bpc) ■ (8) 

Next, this relation is used to derive the flux terms from cell interior to cell 
wall. Thereby, it is assumed that a fraction /3* of all dislocations within a 
distance VcAt to the 6 walls is woven into them, cf. Fig. 4. Here, At is 
an infinitesimal amount of time. That gives a "collecting" volume of V* = 
6vc/^t{d — wy , neglecting the parts which are (Ai)^ small (dark grey regions 
in Fig. 4). Hence, according to Estrin et al. (1998), the loss of dislocation 



^In some applications, the initial values of dislocation densities may be chosen very 
close to zero (Engels et al., 2012). In this case there is a danger of obtaining negative 
dislocation densitites if the original model (6) is used. 
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Figure 4: Cut of the "elementary cell" used to derive evolution equations (6). The shaded 
regions stand for the "collecting area" of dislocations which are woven into the cell walls. 
The dark grey regions can be neglected. 



^ ALjAt (1) 6/3*% 

K {d-wf bd{i-fy/^ ■ ^ ^ 

As can be seen, the flux I/c/K is not vanishing as Pc ^ 0. As a consequence, 
the loss of pc can continue even when there are no longer any dislocations in 
the interior. Such behaviour cannot be expected for real physical systems. It 
is questionable whether Orowan's equation is applicable in case of very low 
dislocation densities since according to (8), ^ 00 as pc ^ 0. However, the 
dislocation velocity cannot exceed a certain hmit due to drag forces. This 
gives rise to a modified Orowan-type equation where the maximum speed is 
bound to a value I'max- A simple approach used in the current study is given 

by 

- = + ^ ■ (10) 

""c 7c/ [bpc> 

In other words, the travelling time of a mobile dislocation is assumed equal 
to the travelling time predicted by Orowan's relation plus a fixed minimum 



length in the cell interior is 
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value. After a few manipulations, the dislocation velocity reads as: 



7c 



1+ 



^^(Pe,7c). (11) 



Here, "^(pcTc) = [1 + 7c/(Vc'^max)]"^ stands for cut-off function since it 
ensures that Vc — )■ fmax for pc — )■ 0. Thus, the singularity is "cut off". 
Introducing this into the previous derivation of (9) we get 

The shear wave velocity represents a reasonable value for the upper bound of 
the mean dislocation velocity Wmax (cf- Johnston and Oilman, 1959; Shehadeh 
et al., 2005). From a mathematical point of view, any function ^(pc; 7c) with 
the property 



0< ;^ 

ope 



< oo (13) 

Pc=0 



could be chosen as an alternative to (11). However, equation (11) offers 
the advantage of a clear interpretation. Using the non-dimensional number 
X = hpcVraaxlic wc find that for x > 99 the difference between Eq. (11) 
and Orowan's equation (8) is less than 1%. Choosing x = 99) 7c = 0.1/s, h 
from Table 1 and fmax = 3100™/s we obtain pc ~ 10^/m^. Thus, assuming 
typical values for dislocation densities ranging from 10^°/m^ in well-annealed 
crystals up to 10^^/m^ for severely plastically deformed metals (Hull and 
Bacon, 2011), Eq. (11) is practically equivalent to Orowan's equation (cf. 
Fig. 5). Another limiting case where the classical Orowan relation may not 
be applicable is given by high-rate plasticity with very large 7 (Krasnikov 
et al., 2011; Austin and McDoweU, 2012). 

Further, in (Mckenzie et al., 2007) there is a discrepancy concerning the 
deformation variable 7 occuring in (6a) and (6b). The plastic shear rate is 
introduced in the model by using Orowan's equation. It is not clear how 
to connect this quantity with an arbitrary macroscopic multiaxial state of 
deformation. In addition, it is not guaranteed that 7 increases monotonically 
with time. Consequently, the shear strain cannot be considered as a loading 
parameter, and the integration with respect to the time cannot be replaced 
by the integration with respect 7. This problem can be avoided by adopting 
the accumulated plastic strain s instead. With the common multiplicative 
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Figure 5: Logarithmic plot of Orowan's equation (8) compared to the modified version 
given by (10). 

split of the deformation gradient F — F^- Fiwe obtain the inelastic (plastic) 
strain rate tensor 

£i = i(£i + |i^) ^ith L,:=F,-Ff^ (14) 

as the symmetric part of the inelastic velocity gradient Li. The accumulated 
plastic strain^ is then defined as 

pT=t I 

s{t):= s{t) dr with s{t) ^ J l\\Di{T)\\ , (15) 

Jt=0 

where ||X|| denotes the Frobenius norm of a second-rank tensor X, cf. def- 
inition (42). In fact, from the simulation results of Mckenzie et al. (2007) 
it is clear that the accumulated strain was used instead of the shear strain. 
Assuming plastic incompressibility and simple shear deformation in the x — z- 
plane offers a relationship between both quantities. With respect to a Carte- 



''Frcquently, s is also referred to as inelastic arc lentgh, true strain or deformation 
degree. 
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sian coordinate system^ one may consider 



lit) 

1 



(14) 



1 



1 



(16) 



Taking the norm and evaluating (15) for the accumulated plastic strain, we 
obtain 

s=|7|/V3. (17) 

Thus, we replace 7 with \/3s. In general, Orowan's relation will be s = ^ vpb 
with the constant ^ depending on the deformation mode. Here, ^ = l/y/S 
need not to be considered since it can be absorbed by model parameters. 
Consequently, the model exactly fulfills Orowan's equation in the case of a 
plastic simple shear. In general, the kinematic relation s — vpb has to be 
considered as additional constitutive assumption. 

The following considerations reveal another weakness of the model. Let 
us assume no production or annihilation of dislocations. In this case, there 
is just a dislocation flux stemming from the first terms in Eq. (5a), (5b). 
Consequently, the total dislocation density (cf. Eq. (3)) should then remain 
constant: 

d 

A = — {/Pw + (1 - f)Pc} 



dt 

f{Pw - Pc 



+ /Pw+(l-/)Pc = 0. 



(18) 



Substituting the flux terms from (5) into (18) we get /(p„ — Pc) = which is 
not fulfilled in general. This contradiction can be resolved by having a closer 
look at the rate of dislocation density. It is obtained by differentiating (2) 
with respect to time. Though the total volume Vt is constant, the fraction 
= fVt of the cell wall and Vc — {1 — f)Vt of the cell interior are time- 
dependent and we get 



/ 



~ y V V ~ V f^^ ' 



(19) 



The density of a dislocation population changes if the length of dislocation 
segments is altered due to generation or annihilation. Density also changes 



^Note that upright indices such as i are labels, wheras italic indices a, 6, . . . stand for 
Cartesian coordinates x, y, z 
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if there is merely a change of volume occupied by the population. The latter 
case implies that only the dislocation length is constant, but existing dis- 
location segments can still move. Thus, the description of cell refinement 
presented by model (6) can be improved. The derivation of (6) only consid- 
eres the change of dislocation density due to and Lc with a fixed /. Thus, 
the terms —f/fpw and +//(! — f)pc were ignored. For a physically consis- 
tent formulation, the changes due to the evolution of / must be taken into 
account, too. This becomes especially important if the wall volume rapidly 
decreases or increases (cf. Section 4). 

Considering all the proposed modifications, a refined dislocation density 
two-phase composite model is obtained: 



W 



ds bdf ' bf 

- ko exp — - — Pw - 7:r^^ v20a) 

dpc . 6a* 



ds M(1-/)V3 VA^c, J 

Compared to (6) it offers the following advantages: 1) The positivity of the 
dislocation densities is guaranteed. 2) The dependence of the model on the 
macroscopic state of deformation is unambiguously defined by Eq. (15). 3) It 
is physically consistent in such a way that without generation or annihilation 
the total dislocation density pt remains constant. 



3. Parameter identification on cell size measurements 

The modified two-phase composite model (20), which has been proposed 
in the previous section contains a number of parameters (cf. Table 1). The 
fundamental physical constants and processing parameters are known. Fur- 
ther, the material parameters a, /o, /oo defining the evolution of the cell wall 
fraction / can be determined from TEM measurments (Mckenzie et al., 2007). 
Hence, seven material constants arranged in the parameter column matrix 

x^[a*,p*,n^,n^,ko,So,K] (21) 
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are left to be identified. Since they can hardly be measured directly, a refined 
strategy for parameter identification has to be used. Basically, the elements 
of X can be found by demanding that the model prediction for a measur- 
able quantity matches best corresponding experimental data. The modified 
model (20) predicts the evolution of the densities pw, Pc of the dislocation 
populations and thus the total density pt and cell size d as function of the 
accumulated plastic strain s. It is well known that dislocation densities are 
extremely hard to measure. Instead of dislocation densities, the cell size 
measurements^ can be adopted. Thus, an error can be defined as follows: 

n = d{x, Si) - d{si) . (22) 

Here, d{x, Si) denotes the model prediction, d{si) stands for the cell size 
measurements at discrete values Si of accumulated plastic strain. Now, the 
optimal set Xopt of parameters within S (ZMJ can be found by minimizing a 
sum of squared errors as target function: 

N 

*(xopt) = min {*(x)} , *(x) = Yl ^'(^) ■ (23) 

i=l 

Here, N stands for the number of measurement data points. Introducing the 
residual r(x) = [ri(x), r2(x), . . . ,rjv(x)]"'", it is possible to write ^(x) = r^r. 
As a general requirement, the experimental data (here d) must provide suf- 
ficient information in order to reliably identify x. The presented model (20) 
contains parameters which control the response under a change in tempera- 
ture T, hydrostatic pressure p or accumulated plastic strain rate s, respec- 
tively. Therefore, experimental data d{si) at varying values of T, p and s 
is required to determine nc,n^, ko, Sq. The ECAP experiments by Mckenzie 
et al. (2007) were conducted under different values of hydrostatic pressure by 
adding a so-called back pressure pb — OMPa and pb — 200MPa. But neither 
temperature nor strain rate have been varied. Still trying to identify the 
parameters nc,n„, /cq, Sq would lead to strong correlations among them and 
most likely poor results. This can be prevented by merging the corresponding 



^Note that the cell size is not to be confused with the grain size. 
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parameters: 



^_=^oexpf-^) fAV'^" , (24) 
V n^kBTJ \soJ 

. (25) 

Effectively, we get one single parameter for each dislocation population con- 
troling its annihilation rate. This reduces the number of parameters and 
suppresses strong correlations between them. The simplified evolution equa- 
tions now read as 

^ = + ) ^(Pc, + ^ J 

- (^A;w + y^^Pw , (26a) 



ds \bK{ 

Since is pressure-dependent, we introduce one such parameter for pe = 
OMPa and for pe = 200MPa. This way, we are able to consider both data 
sets from Mckenzie et al. (2007) in the residual but have only one additional 
parameter. To further reduce the number of parameters, the scaling constant 
K is also taken from Mckenzie et al. (2007) for = OMPa and = 200MPa. 
Finally, a reduced set x = [a* , (3* , k^^, k^^^^, kc] has to be identified. It is 
beneficial to only search within a subset <S C defined by the ranges: 

0<a\/3*<l , A;wo,A;w,oo,^c>0 . (27) 

These restrictions follow immediately from the physical meaning of the con- 
stants. For numerical reasons it is favorable when all parameters show a 
similar order of magnitude. If this is not a priori the case, they can be scaled 
by characteristic factors Wj: 

X = Wx with W = diag{Wi, W2,...) . (28) 
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The factors can be chosen in such a way that the initial set °x: only contains 
values equal to one. Thus, the scaling factors represent the actual initial 
approximation of the parameters. For the present example the following was 
set: 

°x = [1, 1, 1, 1, 1]'^ and W = diag(0.003, 0.006, 3.0, 3.0, 5.0) . (29) 

Consequently, the sum of squared errors is now a function of the dimen- 
sionless^° and scaled parameters: ^ = ^(Wx) = ^(x). To incorporate 
the restrictions from (27) into (23), an interior- reflective Newton algorithm 
(Coleman and Li, 1996) with trust region method is applied to solve the 
minimization problem. A reasonably small target function value was chosen 
as the abort criterion for the procedure. 
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Figure 6: Cell size measurements by Mckenzie et al. (2007) for the aluminum alloy AA 
6016 and simulated curves using Xopt- 



Unfortunately, the given problem is not well posed. Depending on the 
start values (here W) , different sets of parameters leading to equally suitable 



^"Noto that since the parameters from set x have no dimension, no units appear in the 
matrix W. 
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minimizations of the target function were found. Eventually, the parameter 
set Xopt was obtained using the factors from (29): 

[a* - 0.0015, = 0.0017, = 2.7120, A;^^,,,, = 2.4526, h = 4.9641] . (30) 

The fitting result is depicted in Fig 6. Regarding the annihilation parame- 
ter we find /cwo > ^w2oo which is a reasonable result: in the case of higher 
hydrostatic pressure, the annihilation rate is smaller due to the pressure 
dependence of the non-conservative dislocation climb. Note that the annihi- 
lation rate kc in the cell interior is significantly larger than in the cell wall. 
This point will be discussed further in Section 4.3. Even though the number 
of parameters was reduced to five, strong correlations were still found among 
them. Nevertheless, the parameters are physically meaningful as will be seen 
in Section 5. 

4. Extension to non-proportional loading 

4.I. Experimental findings and motivation 

It is experimentally confirmed that cellular dislocation structures may 
become unstable upon stress reversal (Hasegawa et al., 1986). For pure 
aluminum and some Al-alloys, Hasegawa et al. (1975) found that cell walls 
consisting of tangled dislocations are dissolved and become indistinct shortly 
after a change from tension to compression. At the same time, the total 
dislocation density significantly decreases and its distribution becomes more 
uniform. Due to the subsequent monotonic loading, cell patterns emerge 
again. If this stage of deformation is sufficiently long, the total dislocation 
density increases and exceeds its previous value while the cell size decreases 
again. 

It seems promising to extend the modified two-phase composite model to 
capture effects due to load path changes. In order to make the model sensitive 
to such non-proportional deformation, a continuum mechanics technique is 
required. As a practical tool, the finite strain viscoplasticity model from 
Shutov and KreiBig (2008a) is used to evaluate the state of strain and stress. 
In the following, this model is outlined briefly. 

4-2. Finite strain viscoplasticity model 

Let us consider the material model of finite strain viscoplasticity, which 
was proposed by Shutov and KreiBig (2008a). The kinematics of the model 
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are in essence based on the double multiplicative split of the deformation 
gradient 

F = F,-F, , Fi = Fie-Fii. (31) 

The first decomposition can be motivated by the idea of local elastic im- 
loading. The second decomposition was proposed by Lion (2000) in order to 
represent nonlinear kinematic hardening of Armstrong-Frederick type. Ac- 
cording to this assumption, the inelastic part £1 is decomposed into a conser- 
vative part Fie and a dissipative part Fii. Shutov and Kreifiig (2008a) orig- 
inally formulated the material model with respect to physically meaningful 
intermediate configurations. Along with the well-known right Cauchy-Green 
tensor C — ■ F, four tensors of the same type are introduced: 

Ci = ■ Fi , Cn = Fii^ ■ Fii , (32) 

Qe — Ee ' Ee , Qie — Ee ' Ee ■ (33) 

Here, the overset symbols denote the configurations /C and /C, where the 
tensors operate. In addition to the velocity gradient tensor L .— F- F~^, we 
define two inelastic counterparts 

Li = Fi • Fi-i and Lii = Fii • Fii"^ . (34) 

Taking the symmetric parts we obtain the inelastic strain rate and its equiv- 
alent for the microstructure: 

Di=i(Li + LiT) , £ii=^(LiiH-LiiT) . (35) 

Next, in order to capture the nonlinear isotropic hardening, we split the 
accumulated inelastic arc length s into a conservative and a dissipative part 
by s = Sc + Sd- The free energy of the viscoplastic body is then decomposed 
into three independent summands 

'IpiQe, Qie, Se) = 'IpeliQe) + V'kin(^ie) + '0iso(Se) , (36) 

capturing the macroscopic elastic strain energy ■^ei and the energy stored in 
the microstructure related to kinematic (V'km) and isotropic (^iso) hardening. 
Now the weighted Cauchy stress tensor acting on the intermediate configu- 
ration, the backstress tensor and a scalar stress measure are postulated as 

.S^2p^ , X = 2p^ , R = p^. (37) 
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Here p denotes the mass density in the reference configuration. Evaluating 
the Clausius-Duhem inequality, we get 

(& ■ I - ^) B + (Sie • X) ■■ £ii - i?Se > , (38) 

which suggests introducing two "driving forces" of the Mandei stress tensor 
type: 

S:=C^-S-X , f:=Cie-X (39) 

with X = Fic ■ X ■ FiJ. The difference between the Mandei stress Ce'S and 
the back stress X yields the so-called effective stress, which is considered 
as the driving force of (macroscopic) plastic flow. The inelastic flow on the 
microstructural level is assumed to be driven by the back stresses. The 
evolution equations for the state variables are chosen in such a way that 
the model a priori fuUfills the Clausius-Duhem inequality (cf. Shutov and 
KreiBig, 2008a): 

A = Ai|°/|||°|| , Dii = Aix|° , h = Xi^/l'^R. (40) 

Here, the notation Ai stands for the inelastic multiplier, which describes 
the rate of inelastic deformation. It is determined by the definition of the 
overstress^^ in combination with Perzyna's rule: 

.Hirii-4(---). 4, 

where 0o represents the characteristic overstress and ap stands for the initial 
uniaxial flow stress. Assuming the isotropic real- valued functions ip^i, '4'kin 
and ■0iso to be known and imposing initial conditions, the system of con- 
stitutive equations is closed. In the current study, the elastic potential ■^ei 
corresponds to a compressible Nco-Hooke with a bulk modulus k and shear 
modulus n, and ■?/'kin to the Neo-Hookean ansatz with a modulus c/2. By 
chosing ipiso = 1 ' ^1/2, the well-known Voce rule will be obtained. The ma- 
terial parameters are p > 0, x > 0, /9 > 0, 7 G M, ?] > 0, m > 1, ap > 0, 
k > 0, fjL > 0, c > 0. Shutov and Kreifiig (2008a) transformed the constitutive 



^^This notation is used to avoid confusion. Originally, Shutov and Kreifiig (2008a) 
denoted the over stresses by /. 
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equations to the reference configuration in order to simplify the numerical 
treatment. The implementation of the model was discussed by Shutov and 
Kreifiig (2008b, 2010a). 

4-3. Coupling with the refined two-phase composite model 

A common way to measure the intensity of load path changes consists in 
evaluating the angle between former and current loading direction. Towards 
that end, first the scalar product and the Frobenius-norm of two second-rank 
tensors U and V are given as 

U:V:^U-V'^^ tr(£ • V^) , \\U\\ := ^/UTU . (42) 

Generalizing the definition of the angle between two first-rank tensors yields 
a similar invariant for second-rank tensors (cf. Schmitt et al., 1985): 

, . u -.v u 

Wilson and Bate (1994) captured the change of loading direction by evaluat- 
ing (43) with the "old" plastic strain tensor at the time t and the "new" one 
a,t t + At. Yalcinkaya et al. (2009) specified (43) with the inelastic velocity 
gradient instead: 

f old f new 

cos O = cos <(Li°'^ Li""^) = : , (44) 

whereas Viatkina et al. (2007) uses only the symmetric part of Li, i.e. the 
inelastic strain rate tensor: 

A old A new 

COS 6* = COS < (Di°'^ IDr"') = ^ ■ (45) 

ll^^.newjl 

In contrast to (44), the inelastic spin does not enter this definition. In the 
present study, saving old values of Di can be avoided by exploiting the 

benefits of the viscoplastic model. As typical for kinematic hardening of 
Armstrong-Frederick type, the back stress tensor lags behind the inelastic 
strain rate tensor. Thus, X automatically saves the previous loading direc- 
tion and we can postulate the load path change measure as follows: 

cosO ^ cos<{X,Di) ^ X ■■ Di / {\\X\\ \\Di\\) . (46) 
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Monotonic loading is only present when = whereas proportional non- 
monotonic loading is characterized by 6* = or 6* = vr. In the case of general 
non- proportional loading, O can take any value within [0, tt]. 

Remark 1 There is a physical interpretation of definition (46) . The power 
per unit volume of the back stresses with respect to plastic flow is given as 

p = X ■■ Di = \\X\\ \\Di\\cos<{X,Di) (47) 

and arises from the rate of the free energy related to kinematic hardening ( cf. 

(38);.- 

~pi^k^n=-{C^e■X)■■D,, + X■■D,. (48) 

It follows from (40) that the first summand on the right hand side of (48) 
is always negative. This loss corresponds to microstructural processes which 
are dissipative by nature. In constrast to the second summand: if X and Di 
are "parallel", the power according to (47) is positive and increases the free 
energy. If they are opposing, no energy is stored, but released/dissipated. 

Definition (46) is advantageous compared to (45) because it does not lead 
to a single unit impulse immediately after a load path change. Due to the 
evolution of the back stress tensor, then steadily goes back to zero if the 
subsequent loading is monotonic. Changing the load path is thus modelled 
with a certain "inertia" . This is not only numerically favorable but also be- 
lieved to better capture the intrinsic physical behavior of materials exposed 
to load path changes. Motivated by the experimental findings explained in 
Section 4.1 an evolution equation for the cell wall volume fraction is formu- 
lated. Experimental observations show that in the case of load path change, 
the dislocation cell walls become indistinct (Hasegawa et al., 1986). Still, a 
complete dissolution to a uniform dislocation distribution is rarely observed 
(Yalcinkaya et al., 2009). Hence, the process can be considered as an increase 
of wall thickness and its volume fraction. Accordingly, / should rise shortly 
after the load path change, depending on its intensity. Afterwards, if the 
deformation stays monotonic for a sufficiently long time, the wall volume 
fraction / should decrease again. Furthermore, the evolution of / has to 
fulfill some additional requirements: a) it must be continuous, b) it should 
saturate under monotonic loading, c) it should include (7) as a special case 
and d) / must remain within the range (0, 1). An approach which fulfills a) 
- c) a priori is presented in the current study. It is based on the relation: 

^ = {s^{X,DO-a{f-f^)]s , f\t=o = fo. (49) 
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The second term on the right hand side stands for the shrinkage of the wall 
volume fraction / with progressive deformation and would result in Eq. (7). 
The first summand causes a temporary increase of / as a response to load 
path changes. Here, the parameter 6 controls the amplitude and the function 
^{X,Di) e [0, 1] determines the intensity and duration of the response. To 
ensure requirement d), it is sufficient to show that the slope of / is non- 
positive in the vicinity of / = 1: 



ds 
ds 



<max{5^(X,Di)-a(l-/oo)| =5-a(l-/oo) , 
/=i ^ ^ 

< for 5< a(l - /oo) . (50) 

/=i 

Thus, chosing 5 < a{l — foo) keeps / within the permissible range. The 

additional parameter can be obtained by measuring the wall thickness prior 
to and shortly after load path change. The function determining the response 
to load path changes is assumed in the form: 

M{X,D,) = J{e) {\\X\\ /max(||X||)}' . (51) 

Here, J^(6') G [0, 1] is used to measure the intensity of load path change as 
a function of 0{X^jD]). The remaining part is explained as follows: Since 
X saturates in case of constant loading direction, its magnitude indicates 
how long the monotonic phase has lasted. Assuming that the dislocation 
cell pattern becomes more and more distinct under monotonic deformation, 
\\X\\/ max(||X||) e [0, 1] measures indirectly how well defined the structure 
has become. The more distinct it is, the stronger the dissolution and the 
related increase of the cell wall volume will be^^. The back stresses can 
also be related to the microstructure: Internal stresses of the dislocation 
cell structure are considered to be partially responsible for the Bauschinger 
effect, which is taken into account in the viscoplastic model by an evolv- 
ing X. In case of load reversal these internal stresses support dislocation 
re-mobilization and thereby pattern disruption. Here these processes are 
considered as the origin of the energy release described in Remark 1. Since 
the influence of the back stress tensor is thus twofold, it enters (51) quadrati- 
cally. An adequate estimation of the saturation value of the norm of the back 



^^In case of permanently changing loading directions there will be no response since X 
remains close to zero. 
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stresses is the given by max(||X||) 1/x (cf. Shutov and Kreifiig, 2010b). 
^{O) can be regarded as a "response function" for ||X|| max(||X||) and 
should be determined according to experimental findings. 

Intuitionally, the strongest response could be expected at the highest 
value of the load path change measure, i.e. for — tt. Such a complete 
reversal is typical for the Bauschinger test. However, experimental findings 
(Wilson and Bate, 1994; Viatkina et al., 2003) suggest the maximum dis- 
solution of cells at so-called cross tests where = 7r/2, which states a first 
condition for J^{0). Viatkina et al. (2003) and Yalcinkaya et al. (2009) there- 
fore use the term 1 — | cos 6*1 comparable to the intensity function ^{0) in 
their modeling. Two other questions should also be answered to reasonably 
interpolate the proposed function. Do small load path changes {0 0) have 
an effect? Which intensity Jb is to be expected for the Bauschinger-test with 
(9 = TT? Leaving the exact answer open at this stage, we demand from J^{0): 







: < 0crit 


y{0) = < 


■ 


■.0 = n/2 






: 6> = 7r 



J^{n/2) > y{0) V 6> ^ 7r/2 , (53) 

where ©crit defines a critical load path change which has to be exceeded 
to trigger any regular cell dissolution. Due to a lack of experimental data, 
for simplicity we choose J^b = \, ^crit — tt/5 and interpolate J^{0) in a 
piecewise linear way (Fig. 7). 

Additionally to the geometrical shift of the boundary (related to a cor- 
responding movement of dislocations) a reverse dislocation fiux should to 
be incorporated. It is known (Amodeo and Ghoniem, 1988) that until the 
condensation of cell walls to subgrains, the tangled dislocations in the walls 
can still be activated to become mobile again. Motivated by experimental 
findings (Christodoulou et al., 1986) a load path change is assumed to cause 
such a re-mobilization and make a fraction of immigrate to the cell inte- 
rior. There, fewer obstacles hinder dislocation movement and the chance of 
mutual annihilation is higher. Effectively, this flux from wall to interior in 
combination with a higher annihilation rate inside the cells (cf. parameter set 
(30)) leads to a more intense decrease in p^. Since, by definition, cell walls 
are regions of high dislocation density and the magnitude of is at least one 
order higher than p^, the total dislocation density pt decreases shortly after 
a load path change (cf. Eq. 3). Hence, the modeling of non- proportional 
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Figure 7: Different functions J'{0) measuring the intensity of a load path change. 



loading is in agreement with the experimental oberservations by Hasegawa 
et al. (1975). There, it is reported that the total dislocation density decreases 
shortly after load reversal and the distribution of dislocations becomes more 
uniform. The model thus gives an explanation for this temporary loss in pt. 

The derivation of the flux of dislocations from wall to interior is similar to 
the reasoning in Section 2.2. It is assumed that a fraction F* of dislocation 
segments woven in the wall within a distance v^/St from the 6 walls is re- 
mobihzed (Fig. 8). That gives a volume of = &v^/^t{d — wY, neglecting 
the parts which are (Ai)^ small (dark grey regions in Fig. 8). Applying the 
modifled version of Orowan's equation to the re-mobilized fraction of pw, the 
loss of dislocation length in the cell wall is 



(1) 



6r*(i - fy/^.i 



bdf 



(54) 



To let this process set in only after load path changes, the quantity F* has 
to be a function of 0, where O is defined by (46). For consistency, we put 
F* — ^*^{X,Di) with the constant < 7* < 1. Now the model is complete 
and the final set of evolution equations for the three state variables is given 
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Figure 8: Cut of the "elementary cell" of the dislocation pattern. The shaded regions 
stand for the "collecting area" of dislocations from the cell walls which are re-mobilized. 
The dark grey regions can be neglected. 
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(55a) 



(55b) 
(55c) 

(55d) 



In the following this model will be referred to as "load path sensitive two- 
population dislocation cell model". Note that choosing ^ — l^S — Y — ^ 
and provided / — > 0, the original model is included as special case of system 
(55). 
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5. Simulation results 



5.1. Simulation of common Equal Channel Angular Pressing routes 

The goal of the subsequent simulation is to estimate the influence of load 
path changes on the cell size using the example of Equal Channel Angular 
Pressing (ECAP). The most simple operation mode of the ECAP process is 
given by route A (Iwahashi et al., 1998). In this case, the billet is pressed 
through the anglular channel over and over in the same fashion. Neglecting 
any curvature of the channel and heterogenity of the resulting deformation, 
the material flow in the active zone corresponds to a simple shear. Hence, 
for the inelastic deformation gradient and the inelastic strain rate we get: 



1 7^ 
1 
1 



(14) 



*-ab 



1 
2 



1 



1 



(56) 



with the monotonic function 7(t). For the sake of simplicity, the measure 
(45) is used in this section to estimate the load path change. According to 
this definition, the angle between D;"^^ and ^i"^'* remains 0. Clearly, route 
A represents successive monotonic loading. 

Let us now recall the load path which corresponds to route C (Iwahashi 
et al., 1998). After each pass the billet is rotated through 180° about its lon- 
gitudinal axis before it is inserted into the input channel again. Accordingly, 
there is a simple shear in the x — z plane with the shear angle increasing at 
each odd pass and decreasing at each even pass. For the inelastic strain rate 
tensor in the "old" (2N - If^ and "new" 27V*^ ECAP pass we obtain 



old 

ia6 



1 



1 



rSnew 



7 
"2 



1 



1 



(57) 



Effectively, this corresponds to a complete load reversal, in other words a 
Bauschinger-test with — tt: 



cos6> = Di"^'^ ■■ D-r^^ I 1 1 A' 



-1 



(58) 



The loading program which corresponds to route B^. is somewhat more com- 
plex (Iwahashi et al., 1998). After each pass the billet is rotated through 
90° around its longitudinal axis. Consequently, there is a simple shear in 
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the X — z plane in the first pass, whereas in the second pass it is the x — y 
plane. In the third pass the deformation of the first pass is reversed and in 
the fourth pass the same is done for the shear from the second pass. Then, 
the procedure starts over. The coefficients of the inelastic strain rate tensors 
for pass 1 to 4 are summarized in the following table: 
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3 


4 


'■ab 
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Oj 




7 
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r 
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7 
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"0 1 0^ 
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The angles between the corresponding ^i"^"^ and follow as: 

12 23 34 41 

COS O — COS O — COS O — COS = 0. 

Thus, we obtain a sequence of cross tests where O = 7r/2. The loading 
program which corresponds to route E is similar (Barber et al., 2004). There 
is a simple shear both in the x — z and x — y plane. In contrast to route B,, 

the deformation of each pass is immediately reversed in the following pass. 
After four passes, the procedure starts over. The coefficients of the inelastic 
strain rate tensors for pass 1 to 4 are summarized in the following table: 
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[a J 
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1 
OJ 
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The angles between the corresponding D^^^^ and °^^ follow as: 



12 34 23 41 

COS — COS — —1^ COS — COS = 0. 

Hence, we obtain a sequence of alternating Bauschinger {0 = tt) and cross 
tests (0 = n/2). 

Even though the kinematic considerations of different ECAP routes are 
idealized, they are useful for developing awareness of basic load path changes 
of route A, C, and E. Neglecting small elastic deformations such that 
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/ and Fi ^ F, all cases are included in the loading program controlled 
by the deformation gradient 



ab 



1 

1 






1 



(59) 



Here, rigid body rotations of the billet are ignored. Thus, any complex 
ECAP process is reduced to a sequence of shears apphed on different shear 
planes. Assuming additionally a right-angled channel, a final shear strain 
of max{Fxy) = max(F-,;2) = 2 is reached in each step. For simplicity, the 
functions F^zit) and F^yit) are piecewise linear. Now, the load path changes 
can be visualized in the {F^z, -F2.y)-space. For route A and C, F^y ~ and 
the loading history can be visualized in one plane. The load path changes 
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Figure 9: Loading program for ECAP route Be with A'' representing the number of ECAP 
passes. 



caused by route Be and E are illustrated in Fig. 9 and 10, respectively. Here, 
load path changes are clearly visible as kinks in {E^zi Fj;^^) -space. 

All four idealized ECAP routes were used to evaluate the load path sen- 
sitive two-population dislocation cell model described by (55) coupled with 
the macroscopic viscoplastic model from Section 4.2. Since the deformation 
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N 

Figure 10: Loading program for ECAP route E with N representing the number of ECAP 
passes. 

was assumed to be homogeneous, it suffices to solve tlie problem at one single 
material point. For a given local deformation history F{t) according to (59), 
the macro and micro material response in the time interval t G [0, T] was ob- 
tained evaluating the viscoplasticity and two-population-dislocation model 
simultaneously. The material parameters which were used are summarized 
in Tables 2 and 3. 



Tabic 2: Macroscopical material parameters used for numerical simulations 



k/MFa 


/x/MPa 


c/MPa 


7/MPa 


o-p/MPa 


m 


ri/s 


x/MPa"^ 


/3 


73500 


28200 


200 


460 


270 


3.6 


2 • 10^ 


0.005 


5 



The physical constants of the micromodel were taken from Table 1. The 
initial and asymptotic value /o, foo of the cell wall volume fraction and the 
scaling constant K were taken from the measurements in Mckenzie et al. 
(2007). Since the parameter a is not given there, it was estimated in such 
a way that / reaches approximately f^o after 8 ECAP passes for Route A. 
Due to a lack of experimental data, the parameter 6 could only be esti- 
mated exploiting the restriction imposed by (50). For the same reason the 
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re-mobilization parameter 7* was adjusted in such a way as to obtain rea- 
sonable results. All remaining parameters were provided by the identified set 
Xopt from Section 3. Since the simulation was performed assuming constant 
temperature T, hydrostatic pressure p and plastic strain rate s, the merged 
annihilation parameters could be used again. No additional back pressure 
was applied such that — /c^o- 

Tabic 3: Microscopical material parameters used for numerical simulations 



a* 


P* 


7* 


a 


5 






0.0015 


0.0017 


0.0500 


0.3500 


0.0700 


2.7120 


4.9641 


K 


^max/™/ s 


/o 


foo 




6.2000 


3100.0 


0.0700 


0.0500 



The initial values Pw|t=o = 10^^/m^ and Pc|t=o = 10^^/m^ were estimated 
in the following way. The resulting total dislocation density corresponds to 
a typical value for well-annealed undeformed crystals and the difference be- 
tween and pc was supposed to be one order of magnitude. The simulation 
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Figure 11: Semi-logarithmic plot of the evolution of dislocation populations simulating 
ECAP Route A. 

results are presented in Fig. 11 - 18. The obtained evolution of the disloca- 
tion densities for route A serves as reference since it corresponds to monotonic 
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loading. In that case the cell wall volume fraction / decreases monotonically. 
Both dislocation populations and with them the cell size d = K/^^^p^ tend 
to saturate (cf. Fig. 12). Hence, the solution of system (55) behaves in 




Figure 12: Evolution of the cell wall volume fraction / (left) and size d (right) for ECAP 
Route A corresponding to Fig. 11. 

a similar way to the solution of the original model (6). In contrast to the 
numerical results reported by Mckenzie et al. (2007) the dislocation densities 
p^{s) and Pc{s) are monotonic functions of s without any change in the sign 
of curvature for Route A. 





Figure 13: Semi-logarithmic plot of the evolution of dislocation populations simulating 
ECAP Route C. 
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Figure 14: Evolution of the cell wall volume fraction / (left) and size d (right) for ECAP 
Route C corresponding to Fig. 13. 



The load path sensitivity and the resulting non-monotonic course of / 
for routes Be, C and E enriches the dynamic behaviour of system (55). 
Perpetual load path changes lead to periodic behavior (hmit cycle in phase 
space) around a stationary point. As intended by the modeling in Section 
4.3 the amplitudes of the oscillation of the microstructural state variables 
Pwl-s), Pcl-s), /(s) depend on the intensity and duration of the load path 
change. 




Figure 15: Semi-logarithmic plot of the evolution of dislocation populations simulating 
ECAP Route B^. 
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Figure 16: Evolution of the cell wall volume fraction / (left) and size d (right) for ECAP 
Route Be corresponding to Fig. 15. 



Let US confront the numerical results with some more experimental find- 
ings. Viatkina et al. (2003) and Yalcinkaya et al. (2009) actually distinguish 
two scenarios of load path changes with different characteristics. In the case 
of complete load reversal {O = tt, proportional non-monotonic loading), a 
disruption of the cells is reported (Viatkina et al., 2003). There, the wall 
volume fraction / does not change, but the cell size d increases. In the case 
of non-proportional loading with = 7r/2 an increase of / but no significant 
change of d is observed. 
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Figure 17: Semi-logarithmic plot of the evolution of dislocation populations simulating 
ECAP Route E. 



35 



The model proposed in this paper cannot distinguish between these two 
ideahzed cases, but it can describe them by trend. The response for 6* = vr (cf. 
Route C) is less intense than for 6* = 7r/2 (cf. Route Be). Thus, / increases 
less in the case of complete load reversal whereas d still increases. Moreover, 
the model offers the advantage that the right hand side of (55) continuously 
depends on load path changes characterized by well-defined tensor invariants. 




Figure 18: Evolution of the cell wall volume fraction / (left) and size d (right) for ECAP 
Route E corresponding to Fig. 17. 

The curves of the predicted cell size show the following result: the less 
intense the load path changes in the loading history the smaller the achiev- 
able cell size. Because the cell size reduction is a precursor of the subgrain 
structure, this means that Route A would produce cristallites of least size. 
However, due to the choice of the material parameters implemented in this 
study, the simulation results can only serve as an initial qualitative estimation 
of these effects. 

Summarizing this subsection we conclude that load path changes accom- 
panying procedures of severe plastic deformation could essentially influence 
the resulting nanocrystalline microstructure. The obtained evolution of dis- 
location populations differs signficantly depending on whether ECAP route 
A, C, Be or E is apphed. 

6. Conclusions and Discussion 

In this paper, an existing dislocation cell model (cf. Estrin et al., 1998; 
Mckenzie et al., 2007) with two dislocation populations was carefully an- 
alyzed, refined and extended. For the refined model, microstructural pa- 
rameters were identified by minimizing a least squares error function. In a 
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further step, the evolution of the microstructural state variables was made 
sensitive to load path changes. In particular, experimental knowledge about 
the dissolution and disruption of dislocation cells was incorporated. Tensor- 
valued internal variables were utilized to capture load path changes. In order 
to provide these quantities, the microstructural model was coupled with a 
phenomenological model of finite-strain viscoplasticity. Thereby, special em- 
phasis was put on defining a physically interpretable measure of local load 
path changes. The resulting set of equations was hence called "load path sen- 
sitive two-population dislocation cell model". The proposed approach was 
then exemplified by the simulation of various load cases as typical for ECAP 
routes. 

The numerical results can now be summarized. It was found that the 
smallest cell size can be achieved in the case of pure monotonic deformation 
according to ECAP route A. Due to the temporary cell dissolution shortly 
after load path changes, ECAP routes C, Be or E could slowdown the shrink- 
age of the cell size. It is worth noting that even though the transient phases 
of cell dissolution seem to be short related to the entire process, they could 
signficantly affect the final cell size and subgrain size respectively. 

In the present study the proposed "load path sensitive two-population 
dislocation cell model" was considered to be a microscopic model. The cou- 
pling to the macroscopic continuum model of viscoplasticity was performed 
in order to account for an aspect of the evolving microstructure. Both models 
cannot explicitly consider the discrete nature of dislocation motion. Noting 
that only the evolution of dislocation densities is considered reveals some 
limitations of the proposed approach. The collective behavior of dislocations 
is much more complicated under real conditions. Hence, other effects could 
also strongly influence the pattern formation and cell size evolution during 
non-proportional loading. However, the development of a practicable ap- 
proach based on fundamental principles of dislocation kinetics which allows 
for a description of the formation and dissolution of dislocation patterns is 
an open problem. Thus, the current study presents a reasonable idealization 
of the aforementioned phenomena as it is able to model large systems in long 
time frames. 

In order to keep the model practically applicable, simplified phenomeno- 
logical assumptions were made. It should be noted that the numerical pre- 
dictions of the evolution of dislocation densities reflect the real behaviour in 
a qualitative rather than in a quantitative way. Nevertheless, the obtained 
results seem reasonable and might provide some insights into the dislocation 
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cell evolution under temporary cell dissolution. Still, it is an open question 
whether the predicted oscillation of the cell size can be observed experimen- 
tally. Adequate experimental investigations should include the measurement 
of cell size evolution shortly after the load path change, which seems to be 
a challenging task. Furthermore, to reliably identify micromechanical pa- 
rameters, corresponding experiments must be designed in such a way that 
each parameter which appears in (55) takes effect. Ideally, a set of mi- 
crostructural measurements at varying temperatures, hydrostatic pressures 
and plastic strain rates from a procedure of severe plastic deformation is 
required. Furthermore, it seems promising to investigate the impact of non- 
proportional loading on dislocation patterns more in detail. In this way, more 
information about the proposed intensity function J^{0) could be gathered. 

The performed numerical simulations reveal an interesting result: under 
the assumptions that the volume fraction of the cell walls increases and a 
reverse dislocation flux from the wall to the cell interior is initiated shortly 
after the load path change, a drop of the total dislocation density is pre- 
dicted. This agrees well with experimental findings (Hasegawa et al., 1975). 
Moreover, this type of behavior was not explicitly introduced into the model, 
but results naturally from the aforementioned assumptions. 

In the current study, there is only an effect from the macroscopic on the 
microscopic state variables, but no dependence of the macro model on the 
micro model. This feedback can be easily introduced by making the phe- 
nomenological equations (40) - (41) dependent on the dislocation population 
densities or the cell wall volume fraction. Additionally, the temperature in- 
put to the load path sensitive two-population model (55) can be provided 
by the thermoviscoplastic model of Shutov and Ihlemann (2011). With the 
help of the coupled model, the solution of initial boundary value problems 
will be possible. In particular, more realistic simulations of severe plastic 
deformation procedures related to heat generation will become possible. 
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